Army  Research  Laboratory 


Viscoelastic  Analysis  of 
Composite  Fl5rwheels  for  Energy 

Storage 

by  Jerome  T.  Tzeng 


ARL-TR-2610 


November  2001 


Approved  for  public  release;  distribution  is  unlimited. 


20011210  024 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position  unless  so 
designated  by  other  authorized  documents. 

Citation  of  manufacturer's  or  trade  names  does  not 
constitute  an  official  endorsement  or  approval  of  the  use 
thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5069 


ARL-TR-2610  November  2001 


Viscoelastic  Analysis  of 
Composite  Flywheels  for 
Energy  Storage 


Jerome  T.  Tzeng 

Weapons  and  Materials  Research  Directorate,  ARL 


Approved  for  pubKc  release;  distribution  is  unlimited. 


Abstract 


A  viscoelastic  analysis  has  been  developed  to  investigate  stress  relaxation 
and  creep  in  a  multilayered  composite  cylinder  subjected  to  rotation.  The 
analysis  accoimts  for  ply-by-ply  variation  of  material  properties,  fiber 
orientations,  and  density  gradients  through  the  thickness  of  cylinders.  A  closed 
form  solution  based  on  the  corresponding  elastic  problem  is  derived  for  a 
generalized  plane  strain  state  in  a  thick-walled,  mtdtilayered  cylinder.  Laplace 
transform  is  then  applied  to  obtain  the  numerical  solution  of  the  viscoelastic 
problem.  This  report  illustrates  the  derivation  of  the  analytical  model  and 
solution  procedure.  A  numerical  simulation  shows  substantial  creep  and  stress 
relaxation  in  thick-waUed  cylinders  at  an  elevated  temperature.  Viscoelastic 
effects  of  composite  can  result  in  a  drastic  change  of  stress  and  strain  profiles  in  a 
cylinder  over  a  period  of  time,  which  is  critical  in  terms  of  structural  durability 
for  applications  such  as  energy  storage  flywheels. 
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1.  Introduction 


A  composite  flywheel  is  in  general  operated  at  high  tip  speed.  Accordingly, 
centrifugal  forces  resulting  from  the  rotation  of  tiie  rotors  generate  very  high 
tensile  stresses  in  the  radial  and  circumferential  directions.  It  is  very  critical  to 
ensure  that  the  stress  profiles  and  rotor  dimensions  do  not  change  during  the 
service  life  of  the  rotor.  However,  pol)aner  matrix  composites  generally  creep 
over  a  long  period  of  time,  especially  at  an  elevated  temperature  [1].  The 
associated  stress  relaxation  in  the  composite  will  result  in  the  variation  of  stress 
profile  and  dimension,  leading  to  a  potential  failure.  The  objective  of  this 
investigation  is  to  develop  an  analytical  method  to  study  the  viscoelastic 
behavior  of  thick-walled  composite  cylinders  subjected  to  rotation.  The  analysis 
can  be  applied  to  the  design  of  flywheel  machinery  constructed  from  composite 
materials. 

Activities  in  the  research  of  viscoelasticity  have  mainly  concentrated  on  isotropic 
materials,  including  the  studies  in  references  [2-10].  Schapery  [7]  investigated 
the  general  formulation  of  linear  viscoelastic  boimdary  value  problems  of 
composite  materials,  including  the  thermal  viscoelastic  problems  for 
thermorheologically  simple  materials,  and  the  applications  of  the 
correspondence  principle.  A  viscoelastic  analysis  developed  recently  by  the 
author  [11-13]  investigated  thick-walled  laminated  composite  cylinders 
subjected  to  thermal  and  mechanical  loads.  This  investigation  extends  the 
viscoelastic  analysis  to  model  rotating  composite  cylinders. 

In  the  following  research,  the  linear  quasi-static  viscoelastic  behavior  of  a 
thick-laminated  composite  cylinder  is  studied.  The  analysis  accoimts  for  ply-by- 
ply  variation  of  properties  and  fiber  orientations.  The  thick  cylinder  is  assumed 
to  be  in  the  absence  of  thermomechanical  coupling  and  in  the  state  of  generalized 
plane  strain  such  that  aU  the  stress  and  strain  components  are  independent  of  the 
axial  coordinate.  Moreover,  due  to  the  nature  of  axisymmetry,  all  the  stress  and 
strain  components  are  also  independent  of  the  circumferential  coordinate.  The 
mechanical  responses  of  this  thick  composite  cylinder  will,  tiierefore,  only  have 
to  satisfy  the  governing  equation  in  the  radial  direction. 

Invoking  the  Boltzmann  superposition  integral  for  the  complete  spectrum  of 
increments  of  anisotropic  material  constants  with  respect  to  time,  the  viscoelastic 
constitutive  relations  of  the  anisotropic  composite  cylinder  can  be  derived  in 
integral  forms.  Since  the  thick  composite  cylinder  is  maintained  at  a  constant 
elevated  temperature  and  boundary  conditions  are  all  independent  of  time, 
formulations  of  the  linear  thermal  viscoelastic  problem  can  have  forms  identical 
to  those  of  tile  corresponding  linear  elastic  problem  by  employing  the 
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elastic- viscoelastic  correspondence  principle.  In  other  words,  all  of  these  integral 
constitutive  equations  reduce  to  the  algebraic  relations  that  are  very  similar  to 
those  developed  for  elastic  media  when  they  are  Laplace-transformed  by  means 
of  the  rule  for  convolution  integrals.  The  elastic  analysis  can  thus  be  used  to 
derive  the  transformed  viscoelastic  solutions  in  the  time  domain. 


2.  Viscoelastic  Formulation 


The  Boltzmann  superposition  integral  of  stress  Oy  (i,j  =  1,2,3)  and  strain 
eij  (i,  j  =  1, 2, 3)  relation  for  an  isothermal  viscoelastic  problem  is 

=  .  (1) 

OT 

where  C^(t)  is  the  relaxation  modulus  dependent  on  temperature  and  time,  t. 
The  temperature  is  assumed  to  be  constant  in  this  study.  The  Laplace  transform 
of  a  function /ff)  is  defined  as 

7  =  7<s)  =  C  €-'f(t)di.  p) 

where  s  is  the  Laplace  transform  variable.  Appl)dng  equation  (2)  with  the 
convolution  rule  to  equation  (1)  reduces  the  integral  constitutive  equations  to  the 
following  algebraic  relations: 


cry 


C  feu. 


where 


p;  ki  —ki 

C  ..-sr.  . 
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(3) 

(4) 


The  preceding  Laplace  transformed  constitutive  equation  (4)  is  similar  to  elastic 
constitutive  relation  according  to  the  correspondence  principle.  The  transformed 

stiffness  can  be  obtained  from  elastic  moduh,  ,  multiplied  by  Laplacian  "s". 

Consider  a  filament-wound  axisymmetric  thick  composite  cylinder  consisting  of 
N  layers  with  the  axial  coordinate  z,  the  radial  coordinate  r,  and  the 
circumferential  coordinate  9,  as  shown  in  Figure  1.  The  composite  cylinder  has 
the  inner  radius  a,  the  outer  radius  b,  and  the  length  L.  Accordingly,  there  is  a 
corresponding  elastic  problem  with  the  transformed  displacement  components 
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Figure  1.  Coordinate  system  and  stress  components  used  in  a  rotating  cylinder. 

u,  V,  and  w  in  the  axial  direction,  the  circumferential  direction,  and  the  radial 
direction,  respectively,  in  each  layer.  The  axisymmetric  character  of  the  thick 
composite  cylinder,  along  with  the  assumption  of  the  state  of  generalized  plane 
strain,  leads  to  a  simplified  displacement  field,  which  reflects  the  circumferential 

independence  and  only  radial  dependence  of  w , 


u(r,d,z)  =  u(z),  v(r,d,z)  =  v(z),  and  w(r,6.z)  -  w(r).  (5) 

Since  each  layer  of  the  thick  laminated  cylinder  is  cylindrically  monoclinic  in 
reference  to  global  coordinates,  there  is  no  coupling  between  transverse  shears 
and  other  deformations.  Accordingly,  the  vanishing  shear  traction  boundary 
conditions  and  interface  continuity  conditions  generate  zero  out-of-plane  shear 
traction  and  shear  strains  for  each  layer.  Moreover,  owing  to  the  absence  of 
torsional  deformation,  the  transformed  displacement  components  u  and  v 
become 


and 


v  =  0 , 


(6) 

(7) 
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where  the  constant  quantity  £°  has  the  physical  interpretation  of  transformed 

axial  strain  of  a  layer.  In  fact,  >  according  to  the  present  formulation,  also 
represents  the  transformed  axial  strain  of  the  entire  composite  cylinder.  The 
calculation  of  s°  requires  the  knowledge  of  end  boundary  conditions  and  wiU  be 
given  later.  Likewise,  solving  for  w(r)  requires  the  information  of  transformed 
strain  components,  the  constitutive  equations,  as  well  as  the  equilibrium 
equations. 

The  previously  transformed  displacement  field  gives  the  transformed  strain 
components  in  cylindrical  coordinates 


dw(r)  -  w(r)  _  du  -o  _n 

£rr - ’£e0 - £er~  £zr~  £z~^ 

dr  r  dz 


(8) 


The  unabridged  form  of  the  constitutive  equation  (4)  for  each  layer  in  cylindrical 
coordinates  with  the  radial  coordinate  r  normal  to  the  plane  of  symmetry  is 
expressed  as 
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Furthermore,  from  tlie  previous  discussions  it  can  be  shown  that  two  of  tire  three 
equilibrium  equations  are  satisfied  automatically.  The  only  nontrivial 
equilibrium  equation  is  the  one  m  the  radial  direction: 

^_gni^Crrr- (Tee  QO) 

dr  r  s 

Here,  p  is  the  density  of  composite  and  co  is  angular  velocity  of  rotation. 
Substituting  equations  (6),  (7),  and  (8)  into  (9),  the  transformed  stress 
components  and  are  obtained  in  terms  of  the  transformed  radial 

displacement  w.  Incorporating  the  resulting  arr  and  (jee  functions  with 

equation  (10)  gives  a  nonhomogeneous  Euler  differential  equation  of  w  for  a 
layer 
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(11) 
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,  -2~  (Cn-Cn) 

—  +r - 2  w  = 

dr'  dr 


s  C33 


^3 

S  C33 


where 


(12) 


Solving  equation  (11)  for  w  yields  to  the  complete  solution  as  follows: 


Air''  +^2r"^  +  Wp, 


(13) 


where  Ai  2iid  A2  ^^e  coefficients  to  be  determined  from  botmdary  and 
continuity  conditions.  The  particular  solution  can  be  obtained  in  the  following 
expression: 


Wp 


Cl2~Cl3-Oj^  _  ^  j^3 

C33  ”  C22  ^  ~  A,  C33 


(14) 


The  governing  equation  and  the  solution  represents  for  each  single  layer  in  the 
composite  cylinder.  For  a  multilayered  composite  cylinder,  simultaneous 
equations  will  need  to  be  solved  by  appl5dng  proper  boundary  conditions. 
Finally,  it  is  xmderstood  that  the  initial  condition  of  the  original  viscoelastic 
problem  is  a  displacement-free  state  of  rest.  The  boimdary  condition  is  of  free 
traction  and,  hence,  of  free  transformed  traction  on  both  inner  and  outer  circular 
surfaces: 


cr„  =  (7er  =  a^r  ="  0  dXr  =  a.  b. 
On  both  end  surfaces,  stress  resultants  are  zero: 


o^rdr  =  (y„  =  (7.6  =  0  at  z  =  0,L. 


(15) 


(16) 


The  r,  and  To  are  iimer  and  outer  radii,  respectively,  of  an  arbitrary  layer  (the 
layer).  The  continuity  conditions  at  each  interface  between  two  adjacent  layers 
require  continuous  radial  traction  and  continuous  radial  displacement  at  any 
instant.  Thus,  when  written  in  the  transformed  form,  they  become 
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(17) 


—(k)  _  —(k+I)  _ 

^  rr,o  ^ rri 


and 


=  wf\  (18) 

where  k  =  1,...,N  -  1;  and  subscripts  i  and  o  denote  inner  and  outer  surfaces, 
respectively. 

Accordingly,  the  formulation  accounts  for  ply-by-ply  variations  of  material 
properties  and  temperature  change.  The  matrix  form  nmnerical  solution 
procedure  with  parallel  computing  techniques  resolved  the  complexity  and  time- 
consuming  calculation  procedures  in  the  Laplace  transform  of  a  multilayered 
composite  cylinder. 


3.  Numerical  Results 


A  numerical  simulation  is  performed  using  a  thick  composite  cylinder  with  a 
3-in  inner  diameter  and  a  6-in  outer  diameter.  The  cylinder  is  subjected  to 
rotation  of  50,000  rpm.  Two  la5mp  constructions,  an  all-hoop  woimd  [90]3o,  arid  a 
cross-ply  architecture  [(90)4(0)2(90)4]3  are  considered.  This  is  an  assumed  case 
with  ply  thickness  of  0.1  in.  The  cylinders  are  constructed  from  IM7-graphite/ 
8552-epoxy.  Creep  compliance  of  the  8552  epoxy  neat  resin  at  75  °C  is 
characterized  from  DMA  measurements,  as  shown  in  Figure  2.  The  time- 
dependent  characteristic  of  the  resin  is  tihen  used  for  the  creep  compliance  of  the 
composite.  The  creep  compliance  of  the  composite  normalized  with  time  (hour) 
is  listed  as  follows: 


SUt)  =  St  (/r  »  =  7.552,5  X  70-Vpsi , 

Ss(t)=  5°  (r)"“ ,  5"  =  L555¥xio-Vpsi.  (19) 


"t"  is  normalized  time  in  hour.  Here,  the  subscripts  "s  "  and  "tr"  represent  the 
shear  and  transverse  directions  of  composite  materials,  respectively. 

The  subscript  "s"  represents  shear  property.  The  compliance  in  the  fiber 
direction,  5',  =  5.Px7^?'*/psi,  and  Poisson's  ratios,  v/2  =  V/j  =  P-5,  v^3  =  0.5d,  is 
assxuned  to  be  time  independent.  Accordingly,  the  fiber  direction  properties  are 
assumed  to  be  elastic. 


6 


-5.10 


Reference  Tefnperature-ys.O^C 


“S.lSi 


5“b,20-| 

CM 

c 


»-5.25- 

-J 


-5.30-^ 


-5.35- 


4.a59fc:-B  irrj>/lb 


Log  [Time  (hour)  1 


^  DMASuprAs 

Figure  2.  A  DMA  creep  compliance  master  curve  of  8552  epoxy. 


Figures  3-5  show  the  radial  displacement,  radial  stress,  and  hoop  stress 
distribution  in  the  all-hoop  wound  thick-waUed  cylinder  at  three  time  instants, 
respectively.  These  time  instants  noted  in  the  figures  are  initial,  10  years  (10^  hr), 
and  infinite  (lO^o  hr).  Creep  of  composite  materials  over  a  period  of  time  can 
result  in  a  significant  change  in  radial  displacement  and  hoop  strain.  The  radial 
displacement  shown  in  Figure  2  illustrates  significant  creep  over  a  period  of 
time,  which  causes  an  increase  of  the  total  radial  growth.  The  radial  stress 
decreases  significantly  over  a  period  of  time  due  to  stress  relaxation. 
Accordingly,  a  significant  redistribution  is  observed  in  the  hoop  stress  shown  in 
Figure  4;  it  is  critical  from  a  design  point  of  view.  The  gradient  of  hoop  stress 
increases  due  to  viscoelastic  response.  The  hoop  stress  decreases  at  the  inner 
radius  while  it  increases  at  the  outer  radius  over  a  period  of  time.  Accordingly,  a 
safety  margin  has  to  be  imposed  in  the  design  process. 

The  results,  including  the  creep  of  radial  displacement  and  stress  relaxation  of 
radial  and  hoop  stresses,  are  interesting  and  worth  further  examination.  The 
viscoelastic  behavior  of  composite  results  in  more  compliant  radial  and  axial 
stiffriess  in  an  all-hoop  wotmd  cylinder.  However,  the  hoop  stiffness  remains  the 
same  since  it  is  assumed  to  be  elastic.  Figure  6  shows  a  free  body  diagram 
analysis  of  a  section  of  cylinder  that  is  divided  into  two  parts  along  the 
circumference  (not  in  ratio).  Accordingly,  each  free  body  has  to  be  self-balanced 
among  hoop  stress,  centripetal  force,  and  radial  stress.  The  radial  stress,  cr, ,  is  in 
tension  due  to  rotation,  as  shown  in  Figure  4.  The  free  body  at  the  outer  radius 
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Radial  Displacement  (inch) 


Creep  of  Radial  Displacement 


Figure  3.  Creep  of  radial  displacement  in  an  all-hoop  wound  cylinder. 


Relaxation  of  Radial  Stress 


Figure  4.  Relaxation  of  radial  stress  in  an  all-hoop  wound  cylinder. 


Relaxation  of  Hoop  Stress 


Figure  5.  Relaxation  of  hoop  stress  in  an  all-hoop  wound  cylinder. 

must  satisfy  the  force  balance  as  the  following  equation  (conceptually  without 
considering  the  imit): 


2a»0  sin\j/+a,  -m°r’co-.  (20) 

The  radial  stress  decreases  since  the  radial  stiffness  decreases  over  a  period  of 
time.  The  centripetal  force  (left  hand  term)  remains  a  constant  since  the  rotation 
is  the  same.  Accordingly,  the  hoop  stress  ihcreases  at  the  outer  radius  of  the 
cylinder.  On  the  other  hand,  the  free  body  at  the  inner  radius  has  to  satisfy  the 
force  balance  as  follows: 


2  a‘ s  siny/  =  m‘  r'  co-+  a,.  (21) 

The  centripetal  force  remains  a  constant  while  the  radial  stress  decreases. 
Accordingly,  it  is  straightforward  that  the  hoop  stress  will  decrease  at  the  imier 
radius.  Since  both  the  hoop  and  radial  stresses  decrease,  the  radial  displacement 
at  the  inner  radius  also  decreases.  This  analysis  simply  explains  the  physical 
meaning  of  the  stress  and  displacement  variation  due  to  viscoelastic  behavior. 
How'ever,  an  inelastic  analysis  cannot  replace  a  viscoelastic  analysis  that  involves 
time  dependency. 
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Figure  6.  Force  balance  in  the  free  body  diagram  of  a  cylinder. 

Figures  7-9  show  the  radial  displacement,  radial  stress,  and  hoop  stress 
distribution  in  a  thick-walled  cylinder  with  a  layup  construction  of 
[(90)4(0)2(90)4]3  at  three  time  instants,  respectively.  These  three  instants  are 
defined  in  the  previous  section.  The  radial  displacement  shows  characteristics  of 
creep  similar  to  the  all-hoop  wound  cylinder.  However,  the  deformation  profile 
is  not  as  smooth  because  of  the  change  of  fiber  orientation  through  the  thickness. 
The  radial  stress  also  shows  significant  relaxation  over  a  period  of  time.  The 
effect  of  layup  construction  on  the  radial  stress  profile  is  clearly  illustrated  in  the 
"zigzag"  curves.  In  Figure  7,  the  hoop  stress  profile  shows  significant  variation 
over  a  period  of  time.  The  hoop  stress  is  low  in  the  0°  plies  because  of  low 
stiffness  in  the  circumferential  direction.  The  stress  gradient  also  increases 
significantly  in  the  90°  plies.  The  stress  increases  in  the  outer  radius,  while  it 
decreases  in  the  inner  radius. 
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Figure  9.  Relaxation  of  hoop  stress  in  a  cross-ply  wound  cylinder. 


4.  Conclusions 


An  analysis  has  been  dev'eloped  for  the  viscoelastic  behavior  of  laminated 
composite  rotating  cylinders  with  ply-by-ply  variation  of  anisotropic  viscoelastic 
properties.  Stress  relaxation  and  creep  are  properly  calculated  in  the  cylinders 
with  various  layup  constructions.  Creep  and  stress  relaxation  exist  in  the  fiber 
direction  even  though  the  fiber  dominant  properties  are  elastic.  This  mainly 
results  from  the  contribution  of  the  Poisson's  effects.  Relaxation  of  radial  stress 
due  to  the  viscoelastic  transv'erse  and  shear  p)roperties  of  composite  changes  the 
hoop  stress  {profile  as  well  as  radial  displacement. 

The  results  illustrate  that  the  vdscoelastic  behavior  is  critical  for  the  composite 
flywheel  design  to  serve  a  long  lifecycle.  Tlie  stability  of  radial  dimension  is 
critical  for  applications  such  as  rotating  machines  for  power  generations.  The 
rotor  design  has  to  consider  the  significant  change  in  the  hoop  stress  distribution 
and  radial  growth.  This  iiwestigation  provides  a  theoretic  base  to  predict  and 
understand  the  time  dependent  behavior  of  composite  flywheels.  An  analytical 
tool  is  also  developed,  which  can  be  used  for  composite  rotating  machine  design. 
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